Folded Graphene Nanoribbons with Single and Double Closed Edges 
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Graphene nanoribbon folds with single and double closed edges are studied using density func- 
tional theory methods. Van der Waals dispersive interactions are included via semi-empirical pair- 
wise optimized potential. The geometrical phases of the single and double folded ribbons are ob- 
tained. The electronic structure in terms of energy needed for the folding process, van der Waals 
contribution, energy band gaps, and bandstructures are also calculated. The results are interpreted 
in terms of peculiarities of the structures and dispersion interactions. It is shown that significant 
modifications in the electronic structure can be achieved as a result of folding. 



The recent isolation of a single graphene layer has 
generated much scientific interest [U [2]. With its ex- 
traordinary electrical, optical, and mechanical proper- 
ties, this atomically thick two-dimensional material is 
highly expected to be the material of novel technolog- 
ical devices, especially for high-speed electronics [SHS]. 
Graphene nanoribbons (GNRs), which are quasi-one di- 
mensional strips of graphene, have also been synthesized 
[THQ]. GNRs have their own unique characteristics. It 
turns out that the types of edges and the magnitude of 
the width largely determine the metallic or semiconduct- 
ing nature of GNRs. Magnetic states are also found to 
exist in zigzag-terminated nanoribbons p!QUT2 |. 

Graphene and GNRs can be easily folded in the out- 
of-plane direction. Such wrapping changes their func- 
tionalities and properties. Folded GNRs (FGNRs) can 
be viewed as two flat sheets connected by a curved edge. 
They are also referred to as fractional nanotubes because 
the curved edge has a nanotube-like structure. Graphene 
folds have been observed experimentally [13], and they 
can also be achieved by mechanical stimulations [Mj [15] 
or by high temperature annealing [16 . Recently, a vari- 
ety of different forms of graphene pleats have also been 
demonstrated [T7l|18]. The majority of natural graphite 
occurs in the Bernal (AB) stacking sequence, and only a 
small portion appears in the rhombohedral ABC form 
[19 . However, scanning tunneling microscopy studies 
have shown that FGNRs can occur not only in AB^ but 
also in AA or other patterns [I6l [17] . Thus folding can 
offer the possibility to generate other registry dependent 
configurations. 

The existence of interlayer coupling together with the 
edges at the open ends and bending energy at the closed 
edge influence the stability, and electronic and transport 
characteristics in a profound way. Recent studies have 
developed a coarse-grained model showing that the bal- 
ance between the bending and adhesion determines self- 
folding [20 . Nearest-neighbor tight binding reports have 
calculated the conductance of FGNRs showing that the 
electronic structure can be changed significantly by as- 
suming different values of the adhesion between the flat 



portions [2P. Electronic structure calculations have re- 
vealed that FGNRs have permanent electric dipoles [22] . 

In this work we investigate FGNR structures with sin- 
gle and double closed edges. The open sides of the stud- 
ied structures have H saturated armchair ends. The 
folded configurations are found to have a variety of un- 
usual stacking patterns. Our ah initio calculations reveal 
the evolution of the geometries as a function of the width 
of the ribbons. Particular emphasis is placed on the role 
of the van der Waals (vdW) dispersion and its role in 
terms of the stability of the various formations. The elec- 
tronic structure is also obtained showing that single and 
double folds can have much different properties as com- 
pared to those of their planar counterparts. 

Method and Structures. Our calculations are per- 
formed using self-consistent density functional theory 
(DFT) implemented in the Vienna Ab-initio Simulation 
Package (VASP) [23 . The code uses a plane- wave basis 
set and periodic boundary conditions. Here we utilize 
PAW-PBE potentials to account for the core electrons 
[241 [25]. For ah studied systems, (1x1x7) Monkhorst- 
Pack k-grid sampling of the Brillouin zone was taken for 
the self-consistent calculations with an energy cutoff of 
430 eV . All structures are allowed to relax with 10~^ 
eV total energy and 0.02 eV/ A force convergence crite- 
ria. Each supercell is contructed so that the distance 
between the nearest neighbors in each directon is about 
10 i. 

The interlayer coupling is largely determined by the 
long ranged van der Waals (vdW) interaction, which is 
crucial for the formation and stability of the folds. Here it 
is taken into account via the recently developed DFT-D2 
method implemented in VASP [26,. This is a pragmatic 
approach that calculates the vdW interactions originat- 
ing from the fluctuating charge distributions via a semi- 
empirical correction to the conventional Kohn-Sham en- 
ergy. This correction is described with an optimized pair- 
wise force field weighted with a damping function. Full 
parametrization for the entering coefficients for different 
elements was achieved in Refs. [27\ I28j. 

We study armchair edged GNRs, which are nonmag- 
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(c) N=27, >AA-stacking 




(d) A/=28, vAB'-stacking 



FIG. 1: (Color Online) (a) An armchair GNR with H- 
saturated armchair ends for N — 11 Carbon lines. The width 
of the ribbon is W; (b) Front view of the relaxed single FGNRs 
with N = 18, 27, 28, and 35 Carbon lines. Side and top views 
(from left to right in each panel) of the relaxed structures for 
(c) N = 27 - AA stacking; and (d) N = 28 - AB' stacking. 



netic and of semiconducting nature with energy band 
gaps decreasing as the width of the ribbon becomes 
larger. The ends of the ribbons are saturated by H bonds 
in order to avoid chemical self-interaction. Each GNR 
can be represented via the number of C atomic lines TV 
along its axis, as shown in Fig.[l^. A single FGNR is 
prepared by bending the flat ribbon onto itself creating a 
symmetric structure with two parallel graphene portions 
lying on top of each other. A double FGNR is done in 
a similar manner but folding the GNR twice with three 
graphene portions on top of each other. 

Results and Discussion. Each fold is relaxed within 
the above specified criteria. The obtained structures of 
some of the single folded ribbons are shown in Fig.[T]3. We 
find that when 11 < N < 27, the single fold takes a shape 
that looks like a racket without a handle. The distance 
between the furthest points on the curved side increases 
as the width of the ribbon becomes larger with a maxi- 
mum value of 8.7 ^ achieved for N = 27 - Fig.jljD. When 

= 28, the folded structure experiences a geometrical 
phase transition, resulting in a reduced curved region and 
forming two flat GNR-like strips. When A/" > 28, the 
geometrical form becomes racket-like. Interestingly, the 
shape of the curvature is almost unchanged except the 
length of the flat regions becomes larger for larger value 
of A/", as can be seen from Fig.[T]3. The distance between 
the flat regions is similar to the interplane separation in 
graphite. It is 3.2 A for even TV and 3.5 A for odd TV 
number of Carbon lines. 

The type of open edges determines the stacking pat- 
terns of the folded ribbons. It is realized that AA stack- 
ing is geometrically compatible with armchair GNRs, 
while AB stacking is geometrically compatible with 
zigzag GNRs. Fig.[T]3 and Fig.[T]i show the side and top 
views of the studied single folds. It turns out that all rib- 
bons with odd N take AA-stacking patterns. The folded 
ribbons with even N occur in an orthorombic AB' orien- 
tation characterized with a CC bond symmetrically situ- 
ated above the center of each hexagon. We point out that 
AB' stacking is highly unusual in regular graphite. It is 
one of the intermediate phases that may exist if very high 
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FIG. 2: (Color Online) Double FGNRs for: (a) iV = 30 show- 
ing an AAA stacking pattern; (b) N = 35 showing an A A! A!' 
stacking pattern; (c) = 40 showing an AB' A stacking pat- 
tern. 



pressure is applied [29 , and yet here it occurs relatively 
easy by simple folding. 

Fig. [2] shows some of the studied cases of the double 
folded structure. No pronounced geometrical phase tran- 
sition is found here. However, as the GNR increases, the 
flat-like portions become larger, while the curved regions 
experience little change. The distance between the flat 
regions is always 9 ^4, which is almost tripple the dis- 
tance between the flat sides of the single folded ribbons 
with A^ > 28. Considering the top view, the stackings for 
the different ribbons are resolved. It is obtained that they 
fall into categories distinguished by the number of carbon 
lines: for N = 3p {p is integer), the pattern is AAA; for 
A/" = 3p + 1, the pattern is AB'A; and for N = 3p^2, 
it is AA'A". The last geometrical phase is described by 
having two CC bonds symmetrycally positioned in the 
Carbon hexagon (top view). 

The graphene folding is a process that involves a bal- 
ance between the elastic bending and the vdW interac- 
tion. The vdW attraction seeks to make a loop and it 
decreases the overall energy of the structure, while the 
bending tries to resist that tendency and it increases the 
total energy. The balance between these two effects de- 
termines the size and form of the loop as well as the 
flat-like regions. We find that the smallest ribbon for 
which single foldine; can occur is N = 17 {W = 19.7 i), 
while the smallest double folded ribbon has N = 28 Car- 
bon lines {W = 33.3 A). Our calculations show that if 
the vdW potential is not taken into account, the single 
folded N = 17 and double folded A^ = 28 ribbons are 
not stable, indicating the importance of dispersion for 
the stability of the folds. 

Further, we calculate the energy needed to achieve 
folding via the relation 



(1) 



where £^o^i^2 total energies for the unfolded, sin- 

gle folded, and double folded GNRs, respectively. Results 
from the calculations are shown in Fig. [3] It is seen that 
in general, the energy decreases rather smoothly as the 
width becomes larger (larger A^). We compare these re- 
sults with the vdW energy obtained via the D2 Grimme 
approach - Fig.lsb. The vdW interaction is stronger when 
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FIG. 3: (Color Online) (a) Total energy per atom needed to 
create single (Ifold) and double (2fold) FGNRs as a function 
of C atomic lines N; (b) vdW energy per atom for single and 
double FGNRs as a function of C atomic lines N. 
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graphenes or GNRs are parallel and at a distance 3 A. 
Thus Eydw becomes larger as the width of the ribbon 
increases reducing AE. At the same time, the loop is 
smaller, thus less energy is needed for bending. 

For the single FGNR with N = 27, E^^^ is smaller 
in magnitude as compared to the others, which reflects 
the lack of parallel flat-like portions since the geometrical 
phase transition occurs for N = 28. It is interesting to 
note that in practically all of the other cases |£^vdw| is 
in 62 ^ 72 eV range. The oscillatory-like behavior as a 
function of N is attributted to the strong geometrical 
influence originating from the registry dependence of the 
C atoms for the different stackings. This effect has been 
previously realized in other graphitic nanostructures and 
it has also been attributed to the geometrical dependence 
in such dispersive interactions [30 . One notes that for 
single FGNRs with even N {AB')^ £^vdw is always lower 
than that of the closest odd N (AA). The vdW energy 
for the double FGNR with AB' A stacking also appears 
as local minima, however this trend is not pronounced as 
well. 

The electronic structure of the different ribbons is also 
calculated within VASP. The resulting energy bandgaps 
are given in Fig.[4j while the energy bandstructure for 
three GNRs is shown in Fig. [5] All unfolded GNRs are 
found to be semiconductors which is directly related to 
the wave vector quantization along the finite width of the 
ribbon and the shorter H-C bond lengths at the ends. We 
find that the end H-C bonds are about 23% shorter than 
the C-C bonds in the middle of the ribbon (~ 1.42 A). 
This is particularly important for the = 3p + 2 cases, 
which are predicted to be metals within a nearest neigh- 
bor tight binding model [17 . If the tight binding model 
is modified to account for this difference in bond length. 
Eg is non-zero in agreement with DFT calculations. The 
oscillatory and generally decreasing Eg as a function of 
N is in accordance with previously reported ah initio 
simulations of graphene ribbons [31 . The periodicity of 
3 for these band gap oscillations are in the same classes 
as outlined above - N = 3p, 3p + l,3p + 2. These are 
directly related to the tt nature of the ribbon orbitals 
having Fermi wavelength approximately four atomic sites 
along the width of the ribbon. Therefore, every time a C 
line is added, the Fermi wavelength changes in a 3-fold 
periodic pattern. Although the folded structures have 



FIG. 4: (Color Online) Energy band gap Eg for unfolded 
(Ofold), single folded and double folded GNRs as a function 
of C atomic lines N. 



reduced energy gaps as compared to the ones for the un- 
folded ribbons, the oscillatory-like behavior is preserved 
for the single FGNRs. Fig. [4] shows that local maxima 
{N = 3p -\- 1) and minima {3p + 2) in Eg for single folds 
generally coincide with local maxima and minima for the 
unfolds. This behavior, however, is not preserved for the 
double folded structures. This is expected since addi- 
tional stacking patterns occur in the 2fold GNRs. 

The unfolded A" = 28 GNR is a semiconductor with 
relatively large bandgap Eg = 0.461 eV at the F point 
as seen in Fig. [5^. Folding the ribbon once results in 
moving the highest valence and lowest conduction bands 
closer - Fig.[5]3. The energy bands of the double folded 
ribbon around the Fermi level are found to be almost 
touching away from the F point, indicating that the sys- 
tem becomes semimetal. Similar behaviour is found for 
the single folded A^ = 35 ribbon - Fig.[5^. The A^ = 40 
ribbon does not experience a semiconductor-semimetal 
transition due to the folding process. Folding it once and 
twice reduces the gap which is always at F - Fig.|5^,h,i. 
Perhaps the most interesting transition is found for the 
double folded A" = 35 GNR. The system is a metal with 
a band crossing the Fermi level near F - Fig.jsf. 

Our calculations show that the vdW interaction is 
quite important for the magnitude of Eg of the stud- 
ied structures. If the vdW dispersion is not taken into 
account, the separation between the parallel portions is 
larger and consequently the coupling is weaker. For the 
A" = 28 Ifold GNR, for instance, the distance between 
the ends is 3. 356 A without vdW and 3. 156 A with vdW 
interaction, while for the A^ = 40 Ifold GNR, the corre- 
sponding separations are 4.043A and 3.568A. Comparing 
the energy gaps, we find that the vdW coupling reduces 
the band gaps. For example, for N = 39 Eg = 0.138 eV 
(single fold) and Eg = 0.035 eV (double fold) without the 
vdW interaction, while Eg = 0.125 eV (single fold) and 
Eg = 0.0086 eV (double fold) when the vdW interaction 
is taken into account. Similarly, for N = 40 Eg = 0.362 
eV (single fold) and Eg = 0.246 eV (double fold) without 
vdW, while Eg = 0.301 eV (single fold) and Eg = 0.060 
eV (double fold) with vdW interaction. 
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FIG. 5: Energy bandstructures of unfolded, single and double 
FGNRs with N — 28, 35 and 40 Carbon lines, as indicated in 
the top-right corner of each panel. 

Analyzing the electronic structure shows that the en- 
ergy bands around Ep are mainly determined by the 
p orbitals from the edge atoms which are perpendicu- 
lar to GNR axis. These can be considered as edge tt 
states. Upon folding, the interaction between the flat 
parallel-like portions lifts all degeneracies causing the en- 
ergy bands to move closer. If the vdW dispersion is taken 
into account, the interlayer coupling is stronger leading 
to more distorted edge tt states, more prominent band 
splitting and thus smaller gaps. In addition, relatively 
small cr — TT hybridization from the orbitals located on 
the curved regions is found to contribute to the energy 
bands around Ep due to the folding process. The smaller 
distance between the flat-like regions upon inclusion of 
the vdW dispersion in the simulations leads to slightly 
stronger bond changes at the curvature and the a — tt 
hybridization is somewhat increased, although the dom- 
inant effect is still the interlayer interaction. The curved 
regions also have smaller influence on the electronic struc- 
ture of larger ribbons. Their energy bandstructure is 
mainly determined by the parallel regions. The number 
of parallel layers further influences the bandstructure. If 
there are more such layers, the band splitting further 
increases due to the interlayer interaction and further re- 
ducing the band gap. Additional bands may also appear 
around Ep^ which is the case in Fig. 5f, for example. 



Thus, in general, double folded GNRs have smaller Eg. 
For the double folded structures with smaller A^, how- 
ever, there are no significant flat regions since the widths 
of ribbons are relatively small. As a result, for several A^ 
Eg for 2fold is actually larger than Eg for Ifold - Fig.[4j 

These results illustrate how diverse the physical and 
energy band structures can be when effects due curva- 
ture, dispesive interlayer coupling, registry dependence, 
and structure size influence the system simultaneously. 
The stacking patterns and bandstructure characteristics 
functionalities are complex. The shift in registry can re- 
sult in transforming an armchaired GNR from a semicon- 
ductor to another semiconductor, a semimetal, or even a 
metal. Folded ribbons with AA and AB' stackings can 
exist as semiconductors or semimetals depending on the 
number of C lines in each case. Taking into account the 
vdW dispersion is also important since it influences the 
energy band splitting and the resulting energy band gaps. 

Conclusions. Folded graphene nanoribbons with sin- 
gle and double closed edges were considered using first 
principles methods utilizing a DFT-D2 approach. We ob- 
tain geometrical phases of the studied structures in terms 
of various characteristic distances and stacking patterns, 
which are generally not found in open ended layered 
graphenes or GNRs. The energy band gaps and band- 
structures are also obtained showing that the folding pro- 
cess together with relative carbon atomic positions can 
result in semiconductor-semimetal and semiconductor- 
metal transitions. The registry dependence of the dis- 
persive interaction is calculated and it is found that it 
plays a decisive role for the stability and magnitude of 
the energy gaps. 
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